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Abstract 

Photons  emitted  by  a  narrow  laser  b»aa  are  followed  through 
multiple  scattering  events  in  the  ocean  until  registered  by  a  detector 
at  the  source  position.  A  realistic  oc» an  model  Is  used  which  takes 
account  nov.  only  of  molecular  scstteri:  %  (Rayleigh)  and  absorption, 
but  also  scattering  and  absorption  by  the  hydrosols  (Mie) .  The 
single  scattering  function  for  the  hydrosols  is  calculated  from  Mie 
theory  assuming  a  relative  index  of  refraction  of  1.15  and  a  sire 
distribution  with  s  modal  radius  ot  3u .  Targets  with  vsrious  surface 
albedos  (A)  are  introduced  st  various  distances  from  th*>  source.  The 
three  dimensional  path  of  the  photons  la  followed  by  a  Monte  Carlo 
technique.  When  A  i  0.02  the  returned  fl  x  per  c*iit  photon  path  length 
from  the  targsta  is  greater  chan  u-e  background  from  the  laser  beam 
for  any  target  distance.  The  returned  flu*  is  plotted  as  s  function 
of  ths  photon  path  length.  In  practice  tha  detection  distance  is 
limited  by  tha  lowest  flux  which  can  be  detected  and  the  background 
of  natural  light.  Tnhomogeneltlea  in  the  optical  properties  of  the 
ocean  can  alao  be  measured  in  this  way. 


The  authore  are  wiih  the  Department  of  Physics,  Texas  AAH  University, 
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Introduction 

The  distribution  of  the  photons  emitted  in  a  pulse  of  radiation 
from  a  laser  and  subsequently  scattered  by  tint  hydros-.?!  ■  e-  -.—  ’1  as  the 
Rayleigh  scattering  centers  in  the  ocean  provides  important  information 
about  the  properties  and  spatial  distribution  of  these  particles. 

In  addition  there  is  the  practical  problem  of  detecting  an  object 
at  a  given  distance  when  illuminated  by  a  laser  source.  Theoretical 
guidance  in  these  problems  has  not  been  available  since  solutions  of 
the  usual  equations  of  radiative  transfer  are  not  known  when  the  source 
is  a  narrow  beam.  Lldar  measurements  oeneath  the  ocean  have  been  made, 
but  r.o  results  e  been  published  as  yet. 

This  ie  an  ideal  problem  for  the  Monte  Carlo  method,  since  the 
radiation  at  the  detector  can  be  calculated  through  all  orders  of 
multiple  scattering  for  realistic  phase  functions  that  take  account  of 
the  highly  anisotropic  scattering  from  typical  hydrosols.  The  Monte 
Carlo  mathod  used  hare  has  bean  p»evioualy  described  by  us  when  solar 
radiation  la  the  source*  "*  and  when  a  laser  source  ie  used  in  the 
atmosphere.^  The  interaction  of  the  laser  beam  with  water  molecules 
and  hydrosols  is  determined  from  single  scattering  functions,  Tne 
probability  of  scattering  through  a  particular  angle  is  giver  bv  the 
appropriate  single  scattering  function  calculated  ir-i  Mie  theory^  for 
various  polydiaper ae  distributions.  The  scattering  angle  is  chosen 
from  the  appropriate  single  scattering  function  for  each  collision. 

In  the  present  study  the  detector  is  always  located  at  the  position 
of  the  source.  Model  oceans  with  "arlous  amounts  of  turbldltv  are 
considered  as  well  *s  targets  with  a  range  of  surface  albedos. 
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Method 

The  photons  are  emitted  fro*  the  source  uniformly  into  a  core 

with  s  half  angle  of  10  ^  radians.  The  exact  three  dimensional  path 

of  each  photon  la  followed  in  t.,«  ccsputt r  rrrp*r !»«*:-  The  scattering 

angle  Is  chosen  from  the  appropriate  single  scattering  function. 

At  each  collision  a  radiance  is  estimated  at  the  detector.  The  detector 

ie  ajauaed  to  be  located  at  the  position  of  the  source  and  to  have  a 

_■» 

half  width  of  5  *  10  radlana.  At  each  collision  a  choice  is  »ade 
from  an  appropriate  dlatributic.i  function  as  to  whether  the  photon  s 
interacting  with  a  hydroeol  or  a  water  molecule  and  whether  the  photon 
is  scattered  or  ebeotbed. 

A  weight,  initially  unity,  is  associated  with  each  photon.  Colli¬ 
sions  are  forced  o  that  the  photon  never  leaves  the  ocean;  the  weight 
associated  with  the  photon  is  changed  each  time  a  forces  collision  occurs 
so  that  the  correct  result  is  obtained.  A  trajectory  is  also  never 
terminated  because  of  absorption;  instead  the  photon  weight  la  again 
approrr lately  adjusted  to  include  this  possibility.  The  use  of  the 
concept  of  photon  w-ights  is  one  of  several  Important  techniques  which 
are  used  to  reduce  the  variance  of  the  result.  A  photon  trajectory  Is 
terminated  only  when  its  weight  falls  below  a  preassigned  number ,  usually 
taken  as  10 

When  a  photon  undergoes  a  collision  with  a  Rayleigh  scattering 

center,  the  scattering  angle  is  chosen  from  the  appropriate  single 

scattering  function  (-  1  ♦  cos^0).  The  single  scattering  functio"  for 

7  b 

the  hydrosola  was  calculated  exactly  from  Mle  theory  ’  .  The  particulate 
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matter  in  the  water  was  assumed  to  have  -  1.15  and  n2  “  0.001, 
where  n^  and  n£  are  the  real  and  imaginary  parts,  respectively,  of  the 
index  of  refraction  with  respect  to  water.  A  size  distribution  was 
assumed  which  is  constant  for  r  <  lu  and  proportional  to  r^  exp  (-2r) 
for  r  >  ly,  where  r  is  the  radius  of  the  hydrosol.  The  modal  radius 
for  this  particle  size  distribution  is  3y. 

The  single  scattering  function  was  computed  from  the  Mie  theory 
for  particles  with  this  size  distribution  and  index  of  refraction. 

The  result  is  shown  in  Fig.  1  of  Ref.  9.  The  scattering  function  is 
highly  anisotropic  with  strong  forward  scattering.  These  particular 
values  of  the  hydrosol  parameters  were  chosen  after  a  study  of  the 
disparate  values  proposed  in  the  literature*1*  Many  different 
types  of  particles,  both  organic  and  inorganic,  occur  in  ocean  waters 
in  varying  amounts.  The  real  test  of  our  model  lies  in  a  comparison 
of  the  calculated  single  scattering  function  with  experimental  measure¬ 
ments  and  eventually  of  the  calculated  and  measured  radiation  fields. 
Measurements  of  Beardsley*^  of  the  single  scattering  functions  of  nat¬ 
ural  hydrosols  from  ocean  and  coastal  waters  cover  a  range  of  scatter¬ 
ing  angles  from  approximately  10°  to  130°.  Within  that  range  our 
calculated  curve  agrees  well  with  the  measurements  for  ocean  water. 

This  agreement  provides  some  justification  for  the  choice  of  parameters 
in  our  model.  It  also  leads  us  to  hope  that  the  calculated  values  for 
scattering  angles  nearer  0°  and  180°  are  approximately  correct,  regions 
where  unfortunately  there  are  as  yet  no  experimental  measurements. 

Whenever  a  photon  encounters  a  hydrosol  the  new  scattering  angle 
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is  chosen  in  the  Monte  Carlo  calculations  from  a  cumulative  distribu¬ 
tion  function  calculated  from  the  single  scattering  function.  The 
model  is  completed  by  the  specification  of  the  ratios  of  the  atten¬ 
uation  lengths  8  which  determine  the  relative  contributions  of  Rayleigh 
and  Mie  scattering  and  the  fraction  of  radiation  absorbed  by  each  type 
of  scattering  center.  The  following  values  were  assumed;  Bg  r/B,j,  R  » 
0.567,  8g  M/6T  M  »  0.833,  and  for  6T<R/BT  a  range  of  values  from  1.000 
to  0.0714,  where  8g  and  ST  represent  scattering  and  total  attenuation 
lengths,  respectively,  and  R  and  M  represent  the  part  due  to  Rayleigh 
and  Mie  scattering  centers.  These  parameters  should  be  representative 
of  ocean  water  at  X  «*  0.46y  with  varying  degrees  of  turbidity.  The 

ratio  8,_  /B„  determines  the  amount  of  turbidity.  (Note:  in  Ref.  9 

T  j  R  T 

due  to  a  typographical  error  the  value  of  8g  R/BT^R  is  given  incorrectly 
as  0.165;  the  value  used  in  the  calculation  is  0.0165  at  X  =  0.65y.) 

Results 

The  results  were  calculated  first  for  an  ocean  model  with  pure 
Rayleigh  scattering;  this  corresponds  to  an  unusually  clear  ocean  with 
no  appreciable  hydrosol  scattering.  The  returned  flux  as  a  function 
of  the  photon  path  length  is  given  in  Fig.  1.  The  returned  flux  is 
normalized  to  unit  flux  emitted  by  the  source.  In  each  case  the  returned 
flux  is  expressed  as  the  total  flux  returned  for  all  photons  whose 
path  length  lies  within  a  one  meter  path  interval.  Since  the  photon 
path  length  is  the  total  distance  covered  by  the  photon  from  source  to 
detector,  it  is  proportional  to  the  time  delay  between  the  sending  of 
the  pulse  and  its  detection.  The  total  attenuation  length  is  assumed 
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to  be  0.036  m 

The  circles  shown  in  Fig.  1  indicate  the  calculated  returned  flux 
including  all  orders  of  multiple  scattering  when  the  target  is  located 
at  10  m.  The  returned  flux  from  the  target  is  indicated  by  the  circles 
on  a  vertical  line  at  a  photon  path  length  of  20  tn.  The  numbers  beside 
the  circles  indicate  the  albedo  (A)  of  the  target.  In  this  case  when 
A  ■  1,  the  returned  flux  from  the  target  is  several  orders  of  magni¬ 
tude  greater  than  that  returned  by  the  scattering  centers  in  the 
ocean  near  the  target.  Even  when  A  »  0.01,  the  returned  flux  is  signif¬ 
icantly  above  the  oceanic  background. 

Similarly  the  calculated  returned  flux  is  shown  by  other  symbols 

when  the  target  is  at  30  m,  50  m,  and  80  m.  In  each  case  the  returned 

flux  from  the  ocean  for  a  given  photon  path  length  is  nearly  the  same 

regardless  of  the  position  of  the  target.  This  oceanic  background 

decreases  with  photon  path  length  because  of  both  exponential  atten- 
-2 

uation  and  the  r  factor  in  the  intensity.  The  small  fluctuation  in 
the  Monte  Carlo  results  is  evident  upon  comparison  of  these  four  com¬ 
pletely  independent  calculations.  Within  the  assumptions  of  this 
calculation  a  target  could  be  detected  with  equal  ease  at  auy  distance 
provided  only  that  the  low  level  of  the  returned  flu:,  could  be  measured. 
However,  no  other  light  source  than  the  laser  has  been  assumed.  At 
sotte  point  the  returned  flu:;  could  not  be  detected  above  the  ambient 
light  level  in  the  ocean  at  the  detector.  It  is  interesting  that  a 
target  can  be  detected  In  principal  at  any  distance  by  this  method. 

Next  the  influence  of  photon  scattering  from  the  hydrosols  in  a 
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moderately  clear  ocean  was  considered.  The  ratio  8  _/g_  is  chosen 

T  |  R  T 

as  0.5.  Thus  in  half  of  the  collisions  the  photon  interacts  with  a 
hydrosol  and  in  the  other  half  with  a  Rayleigh  scattering  center.  The 
attenuation  length  is  assumed  to  be  0.072  m  The  results  are  given 
in  Fig.  2  for  targets1  at  10,  30,  50,  and  80  m.  The  returned  flux 
decreases  more  rapidly  with  photon  path  length  when  hydrosol  scatter¬ 
ing  is  Included.  Also  the  flux  returned  from  the  target  is  not  only 
smaller,  but  is  closer  to  the  background  from  the  oceanic  scattering. 

The  points  for  a  target  with  A  -  0.01  are  only  slightly  above  the 
background . 

/.  ratio  3t  r/8t  '  1/14  is  assumed  for  the  next  case  in  order  to 

study  the  laser  scattering  in  very  turbid  ocean  or  coastal  water.  The 

results  are  given  in  Fig.  3  for  targets  at  2,  5,  and  10  m.  Since  the 

-1 

attenuation  length  is  now  0.216  m  ,  the  photons  can  not  penetrate  a 
great  distance  in  this  turbid  water.  The  returned  flux  decreases 
much  more  rapidly  with  photon  path  length  compared  to  the  previous  cases. 
Also  the  returned  flux  from  the  target  is  closer  to  the  oceanic  background 
for  a  given  target  albedo.  The  returned  flux  for  a  target  with  A  ■  0.02 
is  only  slightly  above  the  background.  Again  it  is  found  that  the 
returned  flux  from  the  target  can  be  discriminated  from  the  background 
with  equal  ease  no  matter  what  the  target  distance  provided  the  very 
weak  returned  signals  can  be  detected  and  provided  the  background  from 
the  laser  source  is  larger  than  the  background  from  natural  light  sources. 

The  ocean  has  been  assumed  a  homogeneous  medium  in  the  three  cases 
studied  here.  How  do  inhomogeneities  in  the  real  ocean  influence  these 


results'5  Can  cceani"  layers  with  different  optical  properties  be 
distinguished  by  this  method'’  Several  additional  calculations  were 
made  to  study  irhowogeneit ies  in  the  ocean.  In  all  cases  the  target 
was  placed  at  50  s.  For  comparison  the  results  of  the  previous  cal¬ 
culation  With  6^  *  0.5  are  shown  as  squares  in  Fig,  4. 

In  the  first  inhomogeneous  node]  it.  is  assured  that  :T  ’  r.  *  0,5 

i  ,  K  I 

f row  0  <  z  <  2  5  and  free  30  *  z  <  50  a,  and  5_  _/r_  *  0.2  from 

T,R  i 

25  <  z  <  30  n.  The  distance  z  is  measured  from  the  source  and  detector. 
This  model  represents  a  turbid  layer  in  the  ocean  over  the  range 
25  <  z  <  30  ai.  The  results  shown  by  open  circles  shew  that  the  re¬ 
turned  flux  shows  an  appreciable  increase  as  the  photon  path  length 
increases  through  60  a  (the  total  distance  from  source  to  boundary  of 
turbid  layer  back  to  detector).  The  flux  then  decreases  more  rapidly 
than  previously  through  the  turbid  layer.  When  the  photon  oath  length 
is  greater  than  60  a  the  elope  of  the  curve  approaches  its  previous 
value,  but  the  value  of  the  returned  flux  is  now  appreciably  below  the 
corresponding  values  for  the  homogeneous  case  because  of  the  increased 
absorption  by  the  turbid  layer.  This  turbid  layer  noticeably  alters 
the  oceanic  background  and  could  easily  be  detected. 

The  second  inhomogeneous  model  is  the  same  as  the  first  except 
that  the  turbid  layer  is  assumed  from  40  <  z  '  50  ra .  The  returned 
flux  agrees  with  the  homogeneous  model  up  to  a  photon  path  length  of  80  at. 
At  this  point  the  returned  flux  abruptly  increases  and  then  falls  below 
the  homogeneous  model.  The  returned  flux  from  the  target  is  also 
indicated  in  Fig.  4.  In  all  cases  given  here  a  target  with  A  »  0.02 


could  be  distinguished  fros  the  oceanic  backgr.  ..rsd . 

The  third  inhomegenous  model  has  five  different  layers  each  with 
its  own  optical  properties.  The  aim  here  was  to  investigate  the  influence 
of  a  acre  nearly  continuous  variation  in  the  optical  properties  of  the 
ocean.  The  assuaed  properties  of  the  payers  were:  £  / 6,,  *  1/2, 

l  i  K  * 

0  <  z  <  1C  ®;  8,.  / 6  *  1/3,  10  <  z  <  20  m;  B„  D/8_  *  1  / 4 ,  20  <  z  <  30  a; 

i  y  K  T  1  y  K  i 

5  k/6t  1/3,  Y.  <  z  <  40  a:  6_  a/S_  *  1/6,  40  <  z  <  50  m.  At  the 

T .  R  i  T ,  K  i 

first  boundary  (photon  path  length  of  20  as)  the  returned  flux  increases 

at  first  and  the  absolute  value  of  the  slope  of  the  curve  also  increases. 

The  sane  effect  occurs  at  each  boundary,  although  on  a  decreasing 

scale  as  the  relative  change  in  the  properties  of  the  water  is  less  at 

each  boundary.  The  result  is  to  change  the  slope  of  the  returned 

flux  curve  at  each  of  the  boundaries.  A  s&easur eaer.t  of  the  slope  of 

this  curve  provides  information  about  the  local  optical  properties  of 

the  water. 

The  importance  of  multiple  scattering  in  the  calculation  of  the 
returned  flux  depend*  greatly  on  the  parameters  of  the  problem.  In 
general  multiple  scattering  Is  more  important  the  more  turbid  the 
water  end  the  larger  the  photon  path  length.  When  B^  ^/B^  ”  1/2  the 
ratio  of  the  returned  flux  including  multiple  scattering  to  the  single 
scattered  flux  is  about  1.1  at  a  photon  path  length  of  40  m  and  about 
1.5  at  100  m.  When  3^.  g/Sj  "  1/14  the  same  ratio  is  about  1.1  at  a 
photon  path  length  of  4  «,  1.3  at  10  sn,  and  3.0  at  20  m. 
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Conclusion 

When  the  surface  albedo  of  a  target  is  greater  than  about  0.02 
It  can  be  detected  at  any  distance  In  theory  by  a  measuresent  of  the 
returned  flux  fro*  a  laser  as  a  function  of  the  photon  path  length 
(proportional  to  the  delay  time).  In  practice  the  detection  distance 
Is  Halted  by  the  lowest  flux  which  can  be  measured  and  by  the  back¬ 
ground  from  natural  light  sources.  The  slope  of  the  curve  for  the 
returned  flui;  as  a  function  of  the  photon  path  length  varies  with  the 
turbidity  of  the  water.  An  optically  inhomogeneous  occsn  causes 
measurable  changes  in  either  the  value  or  the  slope  of  the  ret  rned 
flux  curve. 

fhla  research  waa  supported  by  the  Office  of  Naval  Research 
through  contract  #  N000  1 4 -68-A-0 308-0002  under  Project  NR  083036. 


References 


1. 

G. 

N . 

Plass 

and 

G. 

W,  Kattawar, 

Appl.  Opt.  7, 

361 

( 1968}  . 

2. 

G, 

N. 

Plaas 

snd 

G. 

W.  Kattawar, 

Appl.  Opt.  7_, 

415 

(19o8)  . 

3. 

G. 

N. 

Plasa 

and 

G. 

W.  Kattawar, 

Appl,  Opt.  7_, 

699 

(1968)  . 

4. 

G. 

W. 

Kattawar  j 

and 

G.  N.  Piaas, 

Appl.  Opt.  7* 

869 

o 

Sh 

CD 

w 

J  . 

G. 

N. 

Plaae 

and 

G. 

W.  Kattawar, 

Appl.  Opt.  % 

2489 

>  (1969^. 

6. 

G. 

N. 

Plaas 

snd 

G. 

U.  Kattawar, 

Appl.  Opt.  (in  press,  1971) 

7.  H.  C.  van  de  Hulst,  Light  Scattering  by  Samll  Particles 
(John  Wiley  &  Sons,  Inc.,  New  York,  1957). 

8.  G.  W.  Kattawar  and  G.  N.  Plasa,  Appl.  Opt.  6,  1377  ,  154.9  (1967). 

9.  G.  N.  Plasa  and  G.  W,  Kattawar,  Appl.  Opt.  8,  455  (1969). 

10.  R.  S.  Fraser  and  W.  H.  Walker,  J.  Opt.  Sc..  Aaer.  58,  1636  (1963). 

11.  Y.  LeGrand,  Ann.  Inst.  Ocean.  19^,  393  (1939). 

12.  E.  0.  Hulbert,  J.  Opt.  Soc.  A«er.  315.  698  (1945), 

13.  W.  R.  G.  Atkina  and  H.  H.  7..  le,  Proc ,  Roy.  Soc.  London  B140, 

321  (1952). 

14.  W.  V.  Burt,  Tcllua  6,  229  (1954). 

G.  F.  Baardaley ,  J.  Opt.  Soc.  A»er.  58,  52  .1968). 


15. 


Figure  1 

The  returned  flux  free  a  laser  source  as  a  function  of  the  photon 

path  length  in  meters .  The  half  angle  of  the  source  and  detector 
-3  -3 

are  10  radians  and  5  x  10  radians  respectively.  Calculated  points 
are  shown  for  the  returned  flux  including  all  orders  of  Multiple 
scattering  for  targets  at  10,  30,  50,  and  80  ■.  The  returned  flux 
from  the  target  is  shown  for  various  assumed  values  of  the  surface 
albedo  (A).  The  returned  flux  is  normalized  pe.  unit  flux  froa  the 
source  and  per  aeter  of  photon  path  length.  Only  Rayleigh  scattering 


la  taken  Into  account. 


RETURNED  FLUX 


Figure  2 


Returned  flux  as  a  function  of  photon  path  length  In  seters.  Both 
Rayleigh  and  hydroscu  scattering  centers  are  ass used .  The  rati^  of 

the  total  Rayleigh  attenuation  length  (8_  )  to  the  total  attenuation 

T 

length  f rn*  both  Rayleigh  centers  and  hydrosols  <V-  8t.«/6t  ■  in- 


See  captlv.<>  to  Fig.  1. 


•  *•00* 


Figure  3 

Returned  flux  ee  a  function  of  photon  path  length  In 
0_  /0_  -  1/1*.  Targets  are  aseuawd  at  2,  5,  and  10 

1  i  X  A 


meters . 
a.  See 


caption  to  Fig.  1. 


*1*.  3 


Figure  4 

Returned  flux  a*  a  function  of  pht  path  lengtn  in  meters.  The 

target  ia  aaeuaed  at  50  c.  Calculated  points  for  a  homogeneous 

ocean  with  8_  /g  "  1/2  are  indicated  by  squares.  Open  circles  art? 

I  i 

the  same  except  for  a  turbid  layer  (6^  R/6T  "  1/5)  fro*  25  *•  r  <  50  e, 
where  s  is  the  distance  from  the  source.  Crosses  are  the  same  model, 
axcept  the  turbid  layer  is  from  40  <  <  50  a.  Solid  circles  indicate 

a  model  with  five  different  layers,  each  with  different  optical  pro¬ 
perties:  e_  _/8_  -  1/2,  0  <  i  <  10  ■;  8.„  _/fi_  -  1/3,  10  <  r  <  20  a, 

I VK  T  i i 

8t  */8t  *  l/4,  20  *  *  «  30  «;  eT  Jt/»T  -  l/s,  30  <  x  <  40  a;  sT  R/eT  -  1/6, 

40  <  i  <  50  a.  Sse  caption  to  Fig.  1. 


